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Abstract 

We study heavy hadron spectroscopy near open bottom thresholds. We employ B and B* mesons 
as effective degrees of freedom near the thresholds, and consider meson exchange potentials between 
them. All possible composite states which can be constructed from the B and B* mesons are studied 
up to the total angular momentum J < 2. We consider, as exotic states, isosinglet states with 
exotic J^^ quantum numbers and isotriplet states. We solve numerically the Schrodinger equation 
with channel-couplings for each state. The masses of twin resonances Zb( 10610) and Zb( 10650) 
recently found by Belle are reproduced. We predict several possible bound and/or resonant states 
in other channels for future experiments. 



PACS numbers: 12.39.Jh, 13.30.Eg, 14.20.-c, 12.39.Hg 



I. INTRODUCTION 



Exotic hadrons are studied extensively in recent hadron physics. There have been 
many analyses which imply that they are multi-quark systems or hadronic molecules. In 
strangeness sector, there are several candidates of exotic hadrons, such as fo(980), ao(980), 
A (1405) and so on. The scalar mesons fo(980) and ao(980) may be regarded as tetra-quark 
systems or KK molecules A (1405) is considered to be generated dynamically by 

KN and vrE |3|]. In charm and bottom sectors, recently many candidates of exotic hadrons 
have been reported in experiments and also actively discussed in theoretical studies 



10|. D,(2317) and Ds(2460) may be tetra-quarks or KD molecules. X(3872), Y(4260), 
Z(4050)^, Z(4250)^, Z(4430)^ and so on are also candidates of exotics states. Especially 
Z(4050)^, Z(4250)^ and Z(4430)^ cannot be simple charmonia (cc) because they are elec- 
trically charged. There are also exotic hadrons in bottom flavors. Yb is the first candi- 
date of exotic bottom hadrons. More recently, Zb(10610)^ and Zb(10650)^ with isospin 
one have been reported by Belle 0, The reported masses and widths of the two res- 
onances are M(Zb(10610)) = 10607.2 ± 2.0 MeV, r(Zb(10610)) = 18.4 ± 2.4 MeV and 
M(Zb(10650)) = 10652.2 ± 1.5 MeV, r(Zb(10650)) = (11.5 ± 2.2) MeV. They also cannot 
be simple bottomonia (bb) because they are electrically charged. 

Well below the thresholds in the heavy quark systems, quarkonia are described by heavy 
quark degrees of freedom, Q and Q (Q = b, c). Above the thresholds, however, it is a non- 
trivial problem whether the resonant states are still explained by the quarkonium picture. 
Clearly, a pair of heavy quark and anti-quark (QQ) are not sufficient effective degrees of 
freedom to form the resonances, because they are affected by the scattering states of the 
two open heavy mesons. Indeed, many resonant states are found around the thresholds in 
experiments. However they do not fit into the ordinary classification scheme of hadrons, such 
as the quark model calculation. Properties for masses, decay widths, branching ratios, and 
so forth, are not predicted by the simple quarkonium picture js]. Therefore it is necessary to 
introduce components other than QQ as effective degrees of freedom around the thresholds. 

Instead of the dynamics of QQ, in the present paper, we study the dynamics described 
by a pair of a pseudoscalar meson P ~ (Qq)spino or a vector meson P* ~ (Qq)spini (q = u, 
d) and their anti-mesons P or P*, which are relevant hadronic degrees of freedom around 
the thresholds. In the following, we introduce the notation P^*) for P or P* for simplicity. 
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We discuss the possible existence of the p(*)p(*) bound and/or resonant states near the 
thresholds. An interesting feature is that the pseudoscalar P meson and the vector P* meson 
become degenerate in mass in the heavy quark limit (Mq — )■ oo). The mass degeneracy 
originates from the suppression of the Pauli term in the magnetic gluon sector in QCD, 



which is the quantity of order (9(1/Mq) with heavy quark mass Mq [11|, |12|. Therefore, 
the effective degrees of freedom at the threshold are given, not only by PP, but also by 
combinations, such as P*P, PP* and P*P*. Because P*^*) includes a heavy anti-quark Q and 
a light quark q, the Lagrangian of P and P* meson systems is give n with respecting the 
heavy quark symmetry (spin symmetry) and chiral symmetry 12|-|21|. 

A new degree of freedom which does not exist in the QQ systems but does only in the 
p{*)p(*) systems is an isospin. Then, there appears one pion exchange potential (OPEP) 
between P'-*^ and P*-*^ mesons at long distances of order l/m,r with pion mass m-,^. What is 
interesting in the OPEP between P*^*) and P^*) is that it causes a mixing between states of 
different angular momentum, such as L and L±2, through its tensor component. Therefore, 
it is expected that the p(*)p(*) systems behave differently from the quarkonium systems. In 
reality in addition to the one pion exchange dominated at long distances, there are multiple 
pion (tttt, tttttt, etc.) exchange, heavy meson (p, w, cr, etc.) exchange at short distances 
as well. With these potentials, we solve the two-body Schrodinger equation with channel- 
couplings and discuss the existence of bound and/or resonant states of p(*)p(*). 

In this paper we study p(*)p(*) systems, with exotic quantum numbers which cannot be 
accessed by quarkonia. The first group is for isosinglet states with / = 0. We recall that 
the possible J^'-^ of quarkonia are J^^ = 0^^ (Vh), 0^^ (Xbo) for J = 0, J (T), J+~ 
(hb), J^^ (Xbi) for odd J > 1, and J , J^^ , J^^ (Xb2) for even J > 2, where examples 
of bottomonia are shown in the parentheses. However, there cannot be J^*" = and 
0^^, J ^ with odd J > 1, and J'^' with even J > 2 in the quarkonia. These quantum 
numbers are called exotic J^*", and it has been discussed that they are the signals for exotics 
including the p(*)p(*) systems and glueballs. The second group is for isospin triplet states 
with / = 1. It is obvious that the quarkonia themselves cannot be isotriplet. To have a 
finite isospin, there must be additional light quark degrees of freedom 221]. In this regard, 
P*^*^ and P'-*^ mesons have isospin half, and therefore the p(*)p(*) composite systems can 
be isospin triplet. We observe that, near the thresholds, the p(*)p(*) systems can access to 
more variety of quantum numbers than the QQ systems. In this paper, we focus on the 
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bottom sector (P = B and P* = B*), because the heavy quark symmetry works better than 
the charm sector. 

In the previous works, Ericson and Karl estimated the OPEP in hadronic molecules 
within strangeness sector and indicated the importance of tensor interaction in this system 



261]. Tornqvist analyzed one pion exchange force between two mesons for many possible 
quantum numbers in 27|, |28||. Inspired by the discovery of X(3872), the hadronic molecular 



model has been developed by many authors 



0, Q, 



29l433l|. For Zb's many works have 



already been done since the Belle's discovery. As candidates of exotic states, molecular 



structure has been studied 



3J-|40|, and also tetraquark structure 4l|-|46|. The existence of 



Zb's has also been investigated in the decays of T(5S') j47H50|. Our study based on the 
molecular picture of p(*)p(*) differs from the previous works in that we completely take into 
account the degeneracy of pseudoscalar meson B and a vector meson B* due to the heavy 
quark symmetry, and fully consider channel-couplings of B*^*-* and In the previous 

publications, the low lying molecular states around Zb's which can be produced from the 
decay of T{5S) were studied systematically and qualitatively (5l|, |52|. Our present work 
covers them also. 

This paper is organized as followings. In section 2, we introduce (i) the tt exchange 
potential and (ii) the npoj potential between B*^*^ and B*^*-' mesons. To obtain the potentials, 
we respect the heavy quark symmetry for the B'^*^B*^*^7r, B'^*^B(*)p and B(*^B'^*)a; vertices. 
In section 3, we classify all the possible states composed by a pair of B(*) and B^*) mesons 
with exotic quantum numbers I'^{J^^) with isospin /, G-parity, total angular momentum 
J, parity P and charge conjugation C . (C in / = 1 is defined only for states of Iz = 0.) 
In section 4, we solve numerically the Schrodinger equations with channel-couplings and 
discuss the bound and/or resonant states of the B^*)B'^*) systems. We employ the hadronic 
molecular picture and only consider the B*^*)B*^*) states. In practice, there are bottomonium 
and light meson states which couple to these states. The effect of these couplings as quantum 
corrections is estimated in section 5. In section 6, we discuss the possible decay modes of 
these states. Section 7 is devoted to summary. 
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II. INTERACTIONS WITH HEAVY QUARK SYMMETRY 



B*^*^ mesons have a heavy anti-quark b and a hght quark q = u, d. The dynamics of 
the B'^*^B*^*^ systems is given by the two symmetries: the heavy quark symmetry for heavy 
quarks and chiral symmetry for hght quarks. These two symmetries provide the vertices of 
TT meson and of vector meson {v = p, u) with open heavy flavor (bottom) mesons P and P* 
(P for B and P* for B*) 

^vHH = -i(3tiHaHbv''{p^)ba + i>^trHaHf,af,^F^^{p)ha, (2) 
where the multiplet field H containing P and P* is defined by 

Ha = ^^ [P:,r - Pa7,] , (3) 



with the four velocity of the heavy mesons ll|. The conjugate field is defined by Ha 



joHl'jo, and the index a denotes up and down flavors. The axial current is given by ~ 



-rdnfc with 




*= ^ . (4) 



where /,r = 135 MeV is the pion decay constant. The coupling constant \g\ = 0.59 for tcPP* 
is determined with reference to the observed decay width F = 96 keV for D* — )■ Dtt {23], 
assuming that the charm quark is sufficiently heavy. The coupling constant g for ttBB* would 
be different from the one for vrDD* because of l/mq corrections with the heavy quark mass 
mq [2^. However the lattice simulation in the heavy quark limit suggests a similar value 



as adopted above 25|, allowing us to use the common value for D and B. The coupling of 
7iP*P*, which is difficult to access from experiments, is also fixed thanks to the heavy quark 
symmetry. Note that the coupling of nPP does not exist due to the parity conservation. 
The coupling constants /3 and A are determined by the radiative decays of D* meson and 
semileptonic decays of B meson with vector meson dominance as /3 = 0.9 and A = 0.56 



GeV ^ by following Ref. 53]. The vector (p and u) meson field is defined by 
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with 



(6) 



and its field tensor by 

Ff,u{p) = df,pu - d^p^ + [p^, p^] , (7) 

where gv = 5.8 is the couphng constant for p — j- tttt decay. 
From Eq. ([1]), we obtain the nPP* and it P*P* vertices 

C^PP* = 2f{PlP:^ + P:lP,)d'^n,,, (8) 

Jtt 

The ttPP* and 'n-P*P* vertices are obtained by changing the sign of the ttPP* and 7[P*P* 
vertices in Eqs. ([8]) and Q. Similarly, from Eq. (|2]) we derive the vPP, vPP* and vP*P* 
vertices {v = p, as 

C,pp = -V2PgvPbPlv ■ pba , (10) 

C^PP. = -2V2\gvv^e^''"'P [pIP^^ - P^P^) d,{po)ba , (11) 

£,p.p, = V2PgvP:P:^v- pha 

+z2V2XgvP:lPUd'{pnba - d'^ip^ba) • (12) 

Due to the G-parity, the signs of vertices for vPP, vPP* and vP*P* are opposite to those 
of vPP, vPP* and vP*P*, respectively, for v = u, while they are the same for v = p. 

It is important that the scatterings p(*)p(*) — ). p(*)p(*) include not only diagonal com- 
ponents PP* — > P*P and P*P* — ). p*p* but also off-diagonal components PP — > P*P* 
and PP* P*P*. The OPEPs for PP* P*P and P*P* P*P* are given from the 



vertices ([8]) and ([9]) in the heavy quark limit as 
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^PlP^^P^P2 



V p* p*_v p* p* 



ri-r2C(r;m^) + S'Ti,T2 T{r;m„) 



and the OPEPs for PP ^ P*P* and PP* ^ P*P* are given as 

V^^P^^p,p. = - {^j;) I [e*,-e*C{r;m^) + S,*,,^^T{r-m^)] t.-t^, 



\'P,P*^P*P* 



(^^2-^ ) ^ ei-f2C{r;m^) + Se*,T2T{r;m^] 



ri-T2. 



(13) 
(14) 

(15) 
(16) 



Here three polarizations are possible for P* as defined by e^"^^ = [^l/\/2,±i/\/2,0) and 
5(0) = (^0, 0, 1), and the spin-one operator T is defined by Ty-^ = ie'^^^e''^ ''^^i'^^- a conven- 
tion, we assign e'^^'> for an incoming vector particle and e^^'>* for an outgoing vector particle. 
Here ti and T2 are isospin operators for P[*^ and -P2*''. We define the tensor operators 



S, 



el,T2 



3(fi-f)(f2-f)-fi-f2, 



(17) 
(18) 
(19) 
(20) 



The p meson exchange potentials are derived by using the same notation of the OPEPs and 
the vertices in Eqs. (fT0l)-(fT2l). 



^PjP2-^PlP2 



^ PIP^^PIP; 



Vp^p*^p*P2 



2m^ 



-C{r;m^)fi-f2, 



(2A(?y)'^ [2£7-£2C(r;m,)-^,.,,,T( 
2 1 r - - 

{2Xgv) - ^2Ti-T2C{r;my)-STj_,T2T{r;m^ 



T1-T2 



+ 



^PiP2^p^p; 
^p^p*^p*p* 



(3gv 

2m^ J 3 
2I 



C(r;m^)fi-f2, 



{2\gv) - [2e*-£2*C(r;m^)-5'£.,£. T(r;m^)] fi-T2, 



{2\gvf - 2e^-f2 C{r; m^)-Sei,T2 T{r;m^) 



TI-T2, 



(21) 

(22) 
(23) 

(24) 

(25) 
(26) 



for V = p. The u exchange potentials are obtained by changing the overall sign from the 
above equations with v = u and by removing the isospin factor ti-T2- 

To estimate the size effect of mesons, we introduce a form factor (A^ — m^)/(A^ + g^) 
in the momentum space at vertices of hPP, hPP* and hP*P* {h = tt, p and u). Here q 
and m/j are momentum and mass of the exchanged meson, and A is the cut-off parameter. 
Then, C(r; rrih) and T(r; rrih) are defined as 



C{r;mh) 



T{r;mh)Si2{r) 



(27r)3 g 2 + ^2 
d^q - 



e*'^-'-F(g;m,), 



{2Tifp + ml 



5i2(g)e*^"^>(g;m,) 



(27) 
(28) 



with Si2ix) = 3{ai-x){a2-x) - ai-(?2, and m?,) = (A2-m^)V(A^ + '? ^)^- The cut-off A 
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is determined from the size of B*^*^ based on the quark model as discussed in Refs. 
There, the cut-off parameter is A = 1070 MeV when the vr exchange potential is employed, 
while A = 1091 MeV when the irpu potential is employed. 

As a brief summary, we emphasize again that, according to the heavy quark symmetry, 
not only the BB* B*B and B*B* B*B* transitions but also the BB B*B* and 
BB* —7- B*B* transitions become important as channel-couplings. In the next section, we 
will see that the latter two transitions supply the strong tensor force, through the channel 
mixing B and B* as well as different angular momentum, such as L and L ± 2. 



III. CLASSIFICATION OF THE b(*)b(*) STATES 

We classify all the possible quantum numbers I^{J^'^') with isospin /, G-parity, total 
angular momentum J, parity P and charge conjugation C for the states which can be 
composed by a pair of B(*) and B^*) mesons. The charge conjugation C is defined for J = 
or J2 = components for 1=1, and is related to the G-parity by G = (— 1)^G. In the 
present discussion, we restrict upper limit of the total angular momentum as J < 2, because 
too higher angular momentum will be disfavored to form bound or resonant states. The 
components in the wave functions for various 

jPC 

are listed in Table HI We use 

the notation "^^^^Lj to denote the total spin S and relative angular momentum L of the 
two body states of B*-*^ and B*^*^ mesons. We note that there are not only BB and B*B* 
components but also BB* ± BB* components. The J^'" = 0^^ state cannot be generated by 
a combination of B^*) and B*^*^ mesons [63] . For / = 0, there are many B(*)B(*) states whose 
quantum number J^'-^ are the same as those of the quarkonia as shown in the third row of 
1 = 0. In the present study, however, we do not consider these states, because we have not 
yet included mixing terms between the quarkonia and the B(*^B'^*^ states. This problem will 
be left as future works. Therefore, for / = 0, we consider only the exotic quantum numbers 
J^*^ = , 1 ^ and 2+~. The states of / = 1 are clearly not accessible by quarkonia. We 
investigate all possible J^*" states listed in Table HI 

From Eqs. (fT3|) -( [T6|) and ( 12T|) -( 126|) . we obtain the potentials with channel-couplings for 
each quantum number I^{J^'^). For each state, the Hamiltonian is given as a sum of the 
kinetic energy and the potential with channel- couplings in a form of a matrix. Breaking of 



8 



TABLE I: Various components of the B(*)b(*) states for several J^*^ ^ 2). The exotic quantum 
numbers which cannot be assigned to bottomonia bb are indicated by ^ . The 0"' state cannot 
be neither bottomonium nor b(*)b(*) states. 





components 


exoticness 


/ = 


I = 1 


0+- 




V 


V 


0++ 


BB(iS'o), B*B*(i5o), B*B*(5Do) 


Xho 


V 


0— 


^ (BB* + B*B) (3Po) 


V 


V 


0-+ 


-)= (BB* - B*B) (3Po), B*B*(3Po) 

V ^ 


Vh 


V 


1+- 


(BB* - B*B) i^Si), (BB* - B*B) {^Di), B*B*{^Si), B*B*(^Di) 


K 


V 


1++ 


(BB* + B*B) i^Si), (BB* + B*B) {^Di), B*B*{^Di) 


Xhi 


V 




BB(iPi), (BB* +B*B) B*B*(iPi), B*B*{^Pi), B*B*(5Fi) 


T 


V 


1-+ 


(BB* - B*B) (3Pi), B*B*(3Pi) 


V 


V 


2+- 


-j^ (BB* - B*B) (3D2), B*B*(3D2) 


V 


V 


2++ 


BB{^D2), (BB* +B*B) (^Ds), B*B*{^D2), B*B*{^S2), B*B*(^D2), B*B*(5G2) 


Xh2 


V 


2-+ 


(BB* - B*B) (3^2), ^ (BB* - B*B) (2^2), B*B*(3p2), B*B*(3F2) 




V 


2— 


^ (BB* +B*B) (3P2), ^ (BB* +B*B) (^^2), B*B*(5P2), B*B*{^F2) 


i^h2 


V 



the heavy quark symmetry is taken into account by mass difference between B and B* mesons 
in the kinetic term. The expHcit forms of the Hamiltonian for each /*^( J^*^) are presented in 
Appendix Rl For example, the J^*" = 1^^ state has four components, 4= (BB* — B*B) (^Si), 
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^ (BB* - B*B) {^Di), B*B*(3^i), B*B*{^Di) and hence it gives a potential in the form of 
4x4 matrix as Eqs. (jA6D, flATfl) and (|A28|) . 



IV. NUMERICAL RESULTS 

To obtain the solutions of the B*^*^B*^*) states, we solve numerically the Schrodinger equa- 
tions which are second-order differential equations with channel- couplings. As numerics, the 
renormalized Numerov method developed in Ref. {sgI is adopted. The resonant states are 
found from the phase shift 5 as a function of the scattering energy E. The resonance posi- 
tion Er is defined by an inflection point of the phase shift S{E) and the resonance width by 



Tr = 2/{d6/dE)E=Er following Ref. 57|. To check consistency of our method with others, 
we also use the complex scaling method (CSM) [58]. We obtain an agreement in results 
between the renormalized Nemerov method and the CSM. 

In Table [Tll we summarize the result of the obtained bound and resonant states, and their 
possible decay modes to quarkonium and light flavor meson. For decay modes, the p meson 
can be either real or virtual depending on the mass of the decaying particle, depending on 
the resonance energy which is either sufficient or not to emit the real state of p or w meson. 
p*{u*) indicates that it is a virtual state in radiative decays assuming the vector meson 
dominance. We show the mass spectrum of these states in FiglH 

Let us see the states of isospin 1 = 1. Interestingly, having the present potential we find 
the twin states in the I^{J^^) = l+(l+~) near the BB* and B*B* thresholds; a bound state 
slightly below the BB* threshold, and a resonant state slightly above the B*B* threshold. 
The binding energy is 8.5 MeV, and the resonance energy and decay width are 50.4 MeV 
and 15.1 MeV, respectively, from the BB* threshold. The twin states are obtained when the 
npu potential is used. We interpret them as the Zb(10610) and Zb(10650) observed in the 
Belle experiment 0, S]. It should be emphasized that the interaction in the present study 
las been determined in the previous works without knowing the experimental data of Zb's 
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Several comments are in order. First, the bound state of lower energy has been obtained 
in the coupled channel method of BB* and B*B* channels. In reality, however, they also 
couple to other lower channels such as nhb, ttT and so on as shown in Table [H Once these 
decay channels are included, the bound state will be a resonant state with a finite width. 
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A qualitative discussion will be given in Section 5. Second, when the vr exchange potential 
is used, only the lower bound state is obtained but the resonant state is not. However, we 
have verified that a small change in the vr exchange potential generates, as well as the bound 
state, the corresponding resonant state also. Therefore, the pion dominance is working for 
the BB* and B*B* systems. (See also the discussion in Appendix[Bl) Third, it would provide 
a direct evidence of these states to be BB* and B*B* molecules if the BB* and B*B* decays 
are observed in experiments. Whether the energies are below or above the thresholds is also 
checked by the observation of these decays. 

In other channels, we further predict the B*^*)B(*^ bound and resonant states. The 
jG^jPC'j _ i-(^o^^) state is a bound state with binding energy 6.5 MeV from the BB 
threshold for the vr exchange potential, while no structure for the npu potential. The ex- 
istence of this state depends on the details of the potential, while the states in the other 
quantum numbers are rather robust. Let us see the results for the latter states from the 
Tcpu) potentials. For 1^(0 ) and 1^(1^^), we find bound states with binding energy 9.8 
MeV and 1.9 MeV from the BB* threshold, respectively. These bound states appear also for 
the vr exchange potential, though the binding energy of the 1^(1++) state becomes larger. 
The 1^(2++) state is a resonant state with the resonance energy 62.7 MeV and the decay 
width 8.4 MeV. The 1~''(1 ) states are twin resonances with the resonance energy 7.1 MeV 
and the decay width 37.4 MeV for the first resonance, and the resonance energy 58.6 MeV 
and the decay width 27.7 MeV for the second. The resonance energies are measured from 
the BB threshold. The l+(2 ) states also form twin resonances with the resonance energy 
2.0 MeV and the decay width 3.9 MeV for the first resonance and the resonance energy 44.1 
MeV and the decay width 2.8 MeV for the second, where the resonance energies have are 
measured from the BB* threshold. 

Next we discuss the result for the states of isospin 7 = 0. In general, the interaction 
in these states are either repulsive or only weakly attractive as compared to the cases of 
1 = 1. The fact that there are less channel-couplings explains less attraction partly. (See 
also Appendix iBl) Because of this, we find only one resonant state with I'~^{J^^) = 0^(1 ^), 
as shown in Fig [1] and in Table UTTI The 0^(1 ^) state is a resonant state with the resonance 
energy 17.8 MeV and the decay width 30.1 MeV for the npu potential. 

In the present study, all the states appear in the threshold regions and therefore are all 
weakly bound or resonant states. The present results are consequences of unique features 
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of the bottom quark sector; the large reduced mass of the B(*)B(*) systems and the strong 
tensor force induced by the mixing of B and B* with small mass splitting. In fact, in the 
charm sector, our model does not predict any bound or resonant states in the region where 
we research numerically. Because the reduced mass is smaller and the mass splitting between 
D and D* is larger. 

TABLE II: Various properties of the B*^*)b(*) bound and resonant states with possible I^{J^^) in 
1 = 1. The energies E can be either pure real for bound states or complex for resonances. The real 
parts are measured from the thresholds as indicated in the second column. The imaginary parts 
are half of the decay widths of the resonances, T/2. In the last two columns, decay channels of a 
quarkonium and a light flavor meson are indicated. Asterisk of p* indicates that the decay occures 
only with a virtual p while subsequently transit to a real photon via vector meson dominance. 





threshold 


E [MeV] 


decay channels 


vr-potential 


irp w-potential 


s-wave 


p-wave 


i+(o+-) 










hb + TT, Xb0,l,2+P 


l-(0++) 


BB 


-6.5 


no 


77b+7r, T+p 


K + P*, Xbl+TT 


i+(o— ) 


BB* 


-9.9 


-9.8 


Xhi+P 


r?b+p, T+TT 


i-(o-+) 


BB* 


no 


no 


hb+/3, Xbo+vr 


T+p 


1+(1+-) 


BB* 


-7.7 


-8.5 
50.4-il5.1/2 


T + TT 


hb + 7T, Xhi+P* 




BB* 


-16.7 


-1.9 


T+p 


hb + P*, XbO,l+TT 


) 


BB 


7.0 - i37.9/2 
58.8 - i30.0/2 


7.1 - i37.4/2 
58.6 - i27.7/2 


hb + TT, Xb0,l,2 + P* 


Vh + P, Tf' + TT 


1-(1-+) 


BB* 


no 


no 


K + P, Xbl+TT 


77b + TT, T + p 


l+(2+-) 


BB* 


no 


no 




hb + TT, Xb0,l,2 + P 


l-(2++) 


BB 


63.5-^8.3/2 


62.7-i8.4/2 


T+p 


hb + P*, Xbi,2 + 7r 


l-(2-+) 


BB* 


no 


no 


hb+P 


T+p 


l+(2— ) 


BB* 


2.0-i4.1/2 
44.2 - Z2.5/2 


2.0 - i3.9/2 
44.1 - Z2.8/2 


Xhi+P 


Vh + P, T + TT 



12 



TABLE III: The b(*)b(*) bound and resonant states with exotic l'~^{J^^) in / = 0. (Same con- 
vention as Table II.) 





threshold 


E [MeV] 


decay channels 


TT-potential 


vr/jw-potential 


s-wave 


p-wave 


o-(o-) 


BB* 


no 


no 




r?b+w, T + r/ 


o+(i-+) 


BB* 


28.6-^91.6/2 


17.8-^30.1/2 






0-(2+-) 


BB* 


no 


no 




K+T], Xb0,l,2+W 



V. EFFECTS OF THE COUPLING TO DECAY CHANNELS 

We have employed the hadronic molecular picture and only considered the B*^*-'B'^*) states 
so far. In reality, however, the B*^*)B*^*) states couple to a bottomonium and a light meson 
state which is predominantly a pion, as Z]Js were discovered in the decay channels of T{nS)T: 



(n = 1, 2, 3) and hb(mP)'7r (m = 1, 2) [J, l5[. In this section, we estimate the effects of such 
channel coupling to the B^*)!?*^*) states. We give a qualitative estimation for the lowest 
B(*)b(*) state in 1^(1+^) corresponding to Zb(10610)^. Similar effects are expected for 
other states. 



To this purpose, we employ the method of Pennington and Wilson [59|. They calculated 
charmonium mass-shifts for including the effect of open and nearby closed channels and 
we apply their calculation procedure for Zb mass-shift. The bare bound state propagator 
i/[s — mg], where mo is the mass of the bare state, is dressed by the contribution of hadron 
loops n(s). Therefore, the full propagator can be written as 

GJs) = J—- = , (29) 

i 



S-ml- En=inn(s) ' 

where s is the square of the momentum carried by the propagator. Ai{s) is the complex 
mass function and the real part of this give the "renormalized" mass. Since the Zb has five 
decay channels, the hadron loops n(s) is a sum of each decay channel n (FigJ2]). Each hadron 
loop n„(s) (Figj2]) is obtained by using the dispersion relation in terms of its imaginary part. 
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B*B*- 
(10650) 



10655 



10649 



10617 



10621 



10622 



Zb(10610) 
gg* __ i 10602 . 

(10604) 10594 

10596 



10606 



10566 

BB 

(10559) 

F(JPC) i+(o--) 1+(1+-) 1-(1^^) 1^(1") 1(2++) l+(2--) 0+(l-+) 



FIG. 1: The b(*)b(*) bound and resonant states with exotic I^{J^^). The dots with error bars 
denote the position of the experimentaly observed Z^'s where M(Zb(10610)) = 10607.2 MeV and 
M(Zb(10650)) = 10652.2 MeV. SoUd hnes are for our predictions for the energies of the bound and 
resonant states when the irp oo potential is employed. Mass values are shown in units of MeV. 

All hadronic channels contribute to its mass at least in principle. Because the dispersion 
integral diverges, we have to subtract the square of mass function A^(.so) at suitable point sq 
from M.{s). Wc shall discuss the choice of Sq shortly. Now, we can write the loop function 
in a once subtracted form as 

An.(., .o) = n.(^) - n„(^o) = /"^"-(/) . (30) 

Jsn [S' - S){S' - Sq) 

Then we arrive at the mass-shift 5M as 

An„(s, So) = M\s) -ml = 5M\s) . (31) 

n=l 

Since an imaginary part of a loop function is proportional to the two-body phase space, we 
take Imll^ in the form for s > s„ as 

Imn„(.) = -gl exp , (32) 
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where Qn is the couphng of Zb to a decay channel n (a bottomonium and a pion), L is the 
orbital angular momentum between a bottomonium and a pion. gem is the magnitude of the 
three momentum of a pion in the center of mass frame and is related to gem by 

qcm=l ,V (33) 



2yi 

In eq f l32p . following [59|], we have introduced the Gaussian-type form factor with a cut- 
off parameter A which is related to the interaction range R. We set A = 600 MeV as a 
typical hadron scale; this value corresponds to -R ~ 0.8 fm by using the relation R ~ \/6/ A. 
Coupling g„ is determined from the partial decay width r„, by r„(s) = —lmlln{s) / y/s . For 
the present rough estimation, we postulate that the decay rates for five final states (T(lS')7r, 
T(2S')7r, T(35')7r, hb{lP)TT, hh{2P)-n) are equal. Then partial decay width for each decay 



channel is set as 3 MeV that is one-fifth of the total decay width 15 MeV 



The subtraction point sq determines the renormalization point where the loop correction 



vanishes. In Ref. j59|, the subtraction point was chosen at the mass of J/ip. Since J /if) 
is a deeply bound state of a cc pair where the charmonium discription works well without 
a DD loop. Now in our situation, there is no such a physical bottomonium like state 
decaying into a pion and a bottomonium. However, as in the case of J/ip we expect that 
the renormalization point of the vanishing loop is located at an energy which is significantly 
below the thresholods of the particle in the loop. We adopt such an energy at ^/sq = 9000 
MeV, 600 MeV below the 7rT(15'), which is similar to the mass difference of J/ip and DD. 

The resulting mass-shift 6M due to each coupling is given in Table IIVI The total mass- 
shift is SM = 2.4 MeV, which is slightly repulsive. This means that the mass of the 3^*^13'^*) 
bound state in 1+(1+^) will be pushed up by the T{nS)TT and hb(mP)7r couplings. Therefore, 
we expect that this state gets closer the BB* threshold, or could even become a resonant 
state. Since the coupling 5'T(i5)7r is the largest due to its low mass, the largest effect is found 
for the coupling of T(lS')7r, where the mass-shift 6M is 6.3 MeV. The coupling of hi,{2P)% 
having P-wave contributes attraction, whose mass-shift 5M is -3.0 MeV. Other coupling 
channels are minor role. 

To summarize this section, we have estimated loop contributions to the mass of the 
B(*)B(*) molecules. We find small repulsive corrections, which still keeps the molecular 
picture unchanged but may change the bound states into resonances, being consistent with 
the experimental observation. 
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r(nS), hJmP) 




FIG. 2: The diagram corresponding to a loop function n„(s) of channel n. 

TABLE IV: Various contributions to loop corrections of channel n, 6M. The total correction is 
shown on the most right column. The first and the second row show the threshold masses and the 
coupling strenghts of channel n. M^h and 6M are given in units of MeV. 





T{IS)ti 


T(25)7r 


T(35)7r 


hb(lP)^ 


hb(2P)7r total 


Mth 


9600 


10163 


10495 


10038 


10399 — 


9n 


1986 


844 


956 


7392 


14179 — 


5M 


6.3 


0.5 


-1.3 


-0.1 


-3.0 2.4 



VI. SEARCH IN DECAYS FROM T(55) 

As twin Zb's were observed from T(5S) decay, T(5S) decay is a useful source to search the 
exotic states around the B*^*)B*^*) energy region. T(5S) can decay to a /"-^(J^*^) = 1"'"(1"'"~) 
state by a single pion emission in s-wave, and a 1^(0 ), 1^(1 ) or l+(2 ) state by a 
single pion emission in p-wave. We recall that the twin Zb's with /"^(J^*^) = 1+(1+") were 
observed in the s-wave channel . In the present study, we further predict the bound 
state in I^{J^'-^) = l+(0 ), and another twin resonant states in I^{J^'-'') = 1^(1 ) and 
1^(2 ) as summarized in Table HTl As for the exotic J^*" states in isosinglet, the resonant 
state in I^{J^'-^) = 0~^{1 can be observed from T{5S) by u emission in p-wave as shown 
in Table [ml 

The radiative decay of T{5S) is also an interesting channel as discussed in Ref. {s^. In 
radiative decay, T(5S) decays to the J^(J^^) = 1"(0++), 1"(1++) and l-(2++) states with a 
photon emission in s-wave. These channels can be also produced in hadronic transitions with 
emission of p meson from higher T-like bottomonim states. In the present study, we predict 
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the bound states in I'^[J^'^) = 1~(0"'"^) and and a resonant states in 1~(2+"'") as 

summarized in Table [Tll 

As a consequence, we will be able to study the B*^*)B*^*) bound and resonant states with 
positive G-parity in a pion emission from T(55') and with negative G-parity in a photon 
emission from T{5S). It will be an interesting subject for experiments to search these states 
in T{5S) decays. 

VII. SUMMARY 

In this paper, we have systematically studied the possibility of the B(*^B'^*) bound and 
resonant states having exotic quantum numbers I^{J^'~^). These states are consisted of at 
least four quarks, because their quantum numbers cannot be assigned by the quarkonium 
picture and hence they are genuinely exotic states. We have constructed the potential of 
the B(*)b(*) states using the effective Lagrangian respecting the heavy quark symmetry. 
Because of the degeneracy in masses of B and B* mesons, the channel mixing, such as BB*- 
B*B, B*B*-B*B*, BB-B*B* and BB*-B*B*, plays an important role to form the B(*)B(*) 
bound and/or resonant states. We have numerically solved the Schrodinger equation with 
the channel-couplings for the B(*)B(*) states with I^{J^^) for J < 2. 

As a result, in / = 1, we have found that the I'^{J^^) = states have a bound 

state with binding energy 8.5 MeV, and a resonant state with the resonance energy 50.4 MeV 
and the decay width 15.1 MeV. We have successfully reproduced the positions of Zb(10610) 
and Zb(10650) observed by Belle. Therefore, the twin resonances of Zb's can be interpreted 
as the B(*)B(*) molecular type states. It should be noted that the BB*-B*B, BB*-B*B* 
and B*B-B*B* mixing effects are important, because many structures disappear without 
the mixing effects. We have obtained the other possible B*^*)B*^*) states in / = 1. We have 
found one bound state in each 1"'"(0 ) and 1~(1''"+), one resonant state in 1~(2+"'") and 
twin resonant states in each 1^(1 ) and 1^(2 ). It is remarkable that another two twin 
resonances can exist in addition to the Zb's. We have also studied the B*^*)B'-*) states in 
/ = and found one resonant state in 0^(1 ^). We have checked the differences between 
the results from the vr exchange potential and those from the npu potential, and found 
that the difference is small. Therefore, the one pion exchange potential dominates as the 
interaction in the B*^*^B*^*^ bound and resonant states. 
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We have estimated the effects of the couphng to decay channels by means of dispersion 
relations. Total mass-shift is 5M = 2.4 MeV, which is slightly repulsive. Therefore, we 
conclude that the molecular picture of B(*)B'^*^ will be a good approximation for the first 
step. More systematic analyses will be left for future works. 

For experimental studies, the T(5S') decay is a useful tool to search the B*^*^B*^*) states. 
T{5S) can decay to the B(*)B(*) states with l+(0"~), and l+(2--) by a single pion 

emission in p-wave and the state with 0^(1 ^) by u emission in p-wave. T{5S) can also 
decay to the B(*)B(*) states with 1~(0++), 1^(1++) and 1^(2++) by radiative decays. In the 
future, various exotic states would be observed around the thresholds from T{5S) decays in 
accelerator facilities such as Belle and also would be searched in the relativistic heavy ion 



collisions in RHIC and LHC \6 



6l| . If these states are fit in our predictions, they will be 



good candidates of the B'^*^B*^*^ molecular states. 
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Appendix A: Hamiltonian 

The hamiltonian is a sum of the kinetic term and potential term as, 

Hjpc = Kjpc + VJpc, (Al) 
for the TT exchage potential only, and 

Hjpc = Kjpc + ^ V}pc, (A2) 

i = TT,p,LO 

for the npu potential. 

The kinetic terms with including the explicit breaking of the heavy quark symmetry by 
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the mass difference uib* — tub are 

Ko++ = diag ( --^Ao, -— + 2AmBB*, --^A2 + 2AmBB* ) , (A3) 
\ 2mBB 2mBB* 2mBB* / 

Ko-+ = diag ( -— ^Ai, -—J^ — Ai + AmBB* ) , (A5) 
\ 2mBB* 2mB*B* / 

Ki+- = diag (--z-r^ — Aq, --prJ^ — A2, --z-J^ Aq + AmBB*, A2 + AmBB* ) , 

\ zmBB* 2mBB* 2mB*B* 2mB*B* / 

(A6) 

=diag ( -— ^Ao,-— ^A2,-— — A2 + AmBB* ) , (A7) 
\ zmBB* 2mBB* zmB*B* / 

Ki— = diag ( — Ao, — Ai + Attibb*, 

V 2mBB 2mBB* 

1 1 1 \ 

Ai + 2AmBB*, -— Ai + 2AmBB*, -^rz A3 + 2AmBB* , (A8) 

2mB*B* 2mB*B* 2mB*B* / 

Ki-+ = diag ( -— ^Ai, -—^ — Ai + AmBB* ) , (A9) 
\ 2mBB* 2mB*B* / 

K2+- = diag ( -— ^A2, -—^ — A2 + AmBB*) , (AlO) 
\ 2mBB* 2mB*B* / 

K2++ = diag ( -— — A2, -—^ — A2 + AmBB*, -TTT^^ ^2 + 2AmBB*, 

\ 2mBB 2mBB* 2mB*B* 

-TT^r^ — Ao + 2 AmBB*, -TT^r^^ — A2 + 2AmBB*,-7rr^ — A4 + 2AmBB* ) , (All) 
^mB*B* zmB*B* ^mB*B* / 

/111 1 \ 

K2-+ = diag Ai, -— A3, -— Ai + AmBB*, -77-: A3 + AmBB* , 

V 2mBB* 2mBB* 2mB*B* 2mB*B* / 

(A12) 

K2— = diag ( — Ai, -773^^ — A3, --—r^ Ai + AmBB*, -77^^^^ A3 + AmBB* ) , 

\ 2mBB* 2mBB* 2mB*B* 2mB*B* / 

(A13) 

where A; = |^ + f |: - with integer / > 0, 1/mBB = I/^b + I/tob, I/ttibb* = 

l/'H^B + 1/'"^B*, l/rnB*B* = l/niB* + 1/mB* and AmBB* = i^^b* — i^^b- 
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The TT exchange potentials for each J^^ states are 



T/TV 
"^0 + + 



^1 



"2+- 



"2- + 



1/1 



Vsvc -Vevx ^ 

a/3Vc 2Vc V2Vt 
— y/eVx V2Vt -Vc + 2Vt ^ 
Vc - 2Vt) , 

Vc + 2Vt -2Vc + 2Vt 
-2Vc + 2Vt Vc + 2Vt 

Vc -V2Vt -2Vc -V2Vt \ 

-V2Vt Vc+Vt -V2Vt -2Vc+Vt 

-2Vc -V2Vt Vc -V2Vt 

-V2Vt -2Vc + Vt -\f2VT Vc+Vt / 

-Vc y2VT VeVr \ 






\/3Vc 



Vc + Vt / 

V3Vc 

-Vc + Vt 

2Vc 

Sa/I^t -^Vc 



(A14) 

(A15) 
(A16) 

(A17) 
(A18) 



-3\/f^T 3^Vr yfvT 

Vc - Vt -2yc - Vt 
-2Vc - Vt Vc - Vt 





-S^Ivt ^ 


3^Vt 


3^Vt 






Vc + |Vt 


— g-Vr 




-Vc + f Vt j 



Vc- Vt 
-2Vc - Vt 




(A19) 

(A20) 
(A21) 



42= Vt 





V 



/ Vc + iVT 



(A22) 



-2^Vt 



3v^ 



Vt 



-2Vc + |Vt 
-^Vt 

/ -Vc - |Vt 
^Vt 
-^Vt 
^Vt 



Vc + |Vt 
-^Vt 



-^vt -vc + fvT; 

-2Vc + iVT 



-Vt \ 



.|Vt 



g-VT -2VC + 

Vc + |Vt -^Vt 

-2Vc + |Vt -^Vt Vc + |Vt / 

3^Vt -^Vt ^Vt \ 

_Vc-|Vt -^Vt ^Vt 

^Vt -Vc-|Vt -|Vt 
-|Vr 



(A23) 



6^Vt 



(A24) 



-Vc + §Vt / 
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The p and ui potentials are 



0++ 



2V^ - 2V^ + V^' 



-4V^. - 2y^ 



2V^ - 2V" + 



( 2Vg + 
V2V^ 



V2V^ 



( 



2Vg -v:^-\- vg' 



V2V^ 



\ s/2v:^ 



-2V3' + Vg" 



-4yg - V 

-v/2y^ 

-2yg + Vr 



V2y^ 
-4yg - 



+ ^c' 

V2yj; 

-\/6y^ 
-v^y,^ 
-2yg - y^ + yg' y 



^yj; 
2yg - yj; + y,^' / 



(A25) 

(A26) 
(A27) 

(A28) 
(A29) 



yv_ 



y;_+ 



"2+- 



y 



2- + 



■ y^j; + y^j; 



V 



2v^yg 







-3 



2V3V^ 













+ Vg' 











6 "^T 


V 6 T 




-2yc" - 



2V^ + V^ + V^' 



2V^ + V^ + V^' 



-4v^ + -y^ 

2V^ +V^ + V^' 
2V^ +V^ + V^' 





2x/3\/(^ 

. IKv" 

V 5 T 
V 7 
35 '^T 

/ 2V^ - + V^' 

5 *'^T 

-4yj - iy- 

3->/6 





-2yg - yj; + v^' 



-syfy^ 



2V3yg 


4Vc^ + yc"' 

§—V'" 



5 "^T 



- V 
5 



7 '^T 



77 '^T 



-2y(5 + 



£ V" -I- T/"' 



35 
12_ 



_^2_yv 



\ aye 
\ 5 



5 "^T 



5 "^T 

ov" — i^y -I- y" 

3-\/6 T/u 
5 T 

-4yg - fyf; 



-4yc- - ^v- 

5 

2VS - iy? + yg' 



_t±-y 

3\/6 T/i; 
~5~ 



-4Vc" - 

3V6 t 
5 ' 



(A30) 

(A31) 
(A32) 



(A33) 



5 T 



(A34) 



2^y^ 2y(^ - |y{; + yg' / 



-2yg + + v',^" 

_3v^T 



3\/6 



y" 



> V" 

6 "^T 
3\/3 T/t) 
5 '^T 

5 "^T 



3 + |yr" + vg' 

3^/2 
5 ''^T 

5 



3V3 
6 

3v^ 



y 

6 T 



2#yn? \ 



-2^0 +I^t"+^c' 



6v/2 

iy 

5 "^T 



-2Vc" 



(A35) 
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where the central and tensor potentials are defined as, 

= [^f^ lc{r;m^)n-T2, (A36) 

= (v^|-) ^T(r;m^)fi-r2, (A37) 

= ~i2\gyY^C{r;m,)n-T,, (A38) 

= {2Xgyf^Cir;m^), (A39) 

V4 = -{2Xgvf^T{r;mp)n-r2, (A40) 

= {2\gyf^T{r;m^), (A41) 

f3gvVi 



Appendix B: Diagonalization of OPEP 



We consider the diagonalization of OPEP ( lA14l) -( lA24p by adopting a stationary approx 



imation for B'-*^ and B^*-* mesons. We regard the B*^*^ and B^*-* mesons as sources of isospin, 
and fix the positions of B^*) and B*^*-* mesons by neglecting the kinetic term. Then the po- 
tentials with channel-couplings as matrices in which have off-diagonal components, turn to 
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be diagonal matrices VJpc as followings, 





= diag 


{-3Vc, -Vc + 4Vt, -Vc - 2Vt) n-r^ 






(Bl) 




= diag 


{-Vc-2VT)n-r2, 






(B2) 




= diag 


{-3Vc,-Vc + Wt)t^-T2, 






(B3) 


Vl+- 


= diag 


i3Vc, 3Vc, -Vc + Wt, -Vc - 2Vt) n 






(B4) 




= diag 


{-Vc + 4W, -Vc - 2Vt, -Vc - 2Vt) 


T1-T2 




(B5) 




= diag 


(3Vc, -Vc + 4Ft, -Vc + 4Vt, -Vc - 


2Vt:, 


-Vc-2V^)ti-T2, 


(B6) 


1 


= diag 


(3Vc,-Vc -2Vt:)ti-T2, 






(B7) 




= diag 


(3Vb,-Vc-2VT)ri-r2, 






(B8) 




= diag 


(3^0, -Vc - 2Vt, -Vc - 2Vt, -Vc - 


2Vt, 


-Vc + 4Vt, -Vc + 4Vt) 




I/2- + 


= diag 


i3Vc, 3Vc, -Vc - 2Vt, -Vc + 4Ft) n 






(BIO) 


V2 — 


= diag 


{-Vc - 2Vt, -Vc - 2Vt, -Vc - 2Vt, 


-Vc 


+ 4V^T)ri-f2, 


(Bll) 



where the central and tensor potentials are defined as, Vc = |C(r;mjr) and Vt = 

{y2j;^^ \T{r] m^), with > and I^t > and Vc < V^. For / = 1 (fi-fs = 1), we see 
that the strongest attractive potential, — (Vc + 2Vt), is contained in the J^^ = 0~^~^, , 
1"' , 1+^, 1 , 1^+, 2+^, 2++, 2^+ and 2 states. In another quantum number, the O""*" 
state in J = 1 has only weakly attractive potential, — 3Vc. Therefore we expect in / = 1 that 
the 0++, 0"", 1"", 1"", 1"+, 2+-, 2++, 2"+ and 2— states may be bound and/or 

resonant states, while the state hardly forms a structure. For 7 = {ti-T2 = —3), the 
strongest attractive potential, —3 ■ 3Vc, is contained in 1 and 2+~ states. The potential 
in the state is repulsive. Therefore there may be a bound and /or resonant states in 
1 ^ and 2+~, and no structure in in J = 0. 

Although the static approximation may be a crude approximation, this is a useful method 
at qualitative level. For example, let us study two nucleon (NN) systems, we analyze the 
deuteron (/ = 0, = 1+) in which the NN potential is given by 2 x 2 matrix with ^5*1 and 
^Di states. [6^] The OPEP for deuteron is given as 



T/NN 



2^/2V^^ 



NN 



2V2V^'' V^"" 



2V^'' 



ri-T2, 



(B12) 



with V^"" 



ffjrNN 
2m ff 



|C(r; m^) and 



NN 



37rNN 

2mjv 



^T(r;m^), for a vrNN vertex constant 
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(^TrNN and a nucleon mass mN- We diagonalize Eq. (1B12P and obtain the eigenvalues in the 
diagonal potential, 

V^^ = diag {V^'' - 4lf ^ V^"" + 2\/r) n-r2- (B13) 

Because > Vq^ , the first eigenvalue gives a repulsion in / = 0, while the second 

eigenvalue gives an attraction in / = 0. The eigenvector of the ^5*1 and ^Di components for 



the second eigenvalue is (a/2/3, a/1/3). This means that the Z^-wave probability is about 33 
percents. In reality, we have to introduce a kinetic term which disfavors the higher angular 
momentum (D-wave), and hence the D-wave probability in deuteron becomes a few (~ 5) 
percents. The stationary approximation will be applied to B*^*^B*^*) systems with better 
accuracy, because the B meson mass is 5.6 times larger than the nucleon mass. We note 
that the OPEP as a long range force are better approximation for larger angular momentum, 
because the sizes of systems are extended. However, the static approximation may become 
worse for larger angular momentum. Therefore, it is necessary for quantitative analysis to 
study numerically the solutions of the Schrodinger equations with the kinetic terms and the 
potentials with channel-coupling as discussed in the text. 
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